Determination of the transverse velocity of Q2237+0305 
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Abstract. The largest systematic uncertainty present in the analysis 
of gravitationally microlensed quasar light curves is that of the galactic 
transverse velocity. We describe a method for determining the trans- 
verse velocity as well as its application to published monitoring data of 
Q2237+0305 (Irwin et al. 1989; Corrigan et al. 1991; 0stensen et al 
1995). We find that while this data displays strong evidence for mi- 
crolensing, it limits the transverse velocity to be < 500A;msec~^ (for 
microlens masses of O.SMq). In addition, by combining these results 
with a new method for modelling the microlensing contribution of stellar 
proper motions, we find that the masses of microlenses in Q2237-I-0305 
are consistent with those of low mass stars. 



1. Introduction 

Q2237-I-0305 is a multiply imaged, gravitationally microlensed quasar. In addi- 
tion to the difficulties presented by the relatively rare occurrence of high magnifi- 
cation events (HME), the analysis of published monitoring data of Q2237-I-0305 
through observation of HME characteristics and rates is hampered by the sparse 
sampling rate. Previous light-curve peaks are under sampled and may have been 
missed altogether. However, comparatively speaking the monitoring data pro- 
vides a good record of the microlensing rate described by the longer timescale 
low level fiuctuations. 



2. Determination of effective galactic transverse velocity 

The observed microlensing rate is produced by the combination of microlens 
proper motions and a galactic transverse velocity. We define the effective trans- 
verse velocity as being that which in combination with a static microlensing 
model, produces a microlensing rate equal to that of the observed light curve 
(Wyithe, Webster & Turner 1999b, c). The effective transverse velocity is there- 
fore larger than the physical transverse velocity. 

We assume that the macro models (eg. Schmidt, Webster &: Lewis 1998) 
for the lensing galaxy in the Q2237-I-0305 system correctly describe the shear 
and optical depth parameters at the position of each image. We assume that the 
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Figure 1. The cumulative histogram for the monitoring data (la- 
belled), and histograms representing the mean (dark line) and ztlu, 
±20" levels (light lines) of fluctuation due to noise only. 

source is small with respect to the microlens Einstein radius (eg Wambsganss 
Paczynski &; Schneider 1990), and find that (in this regime) limits on effective 
transverse velocity and microlens mass are approximately independent of the 
source size assumed. 

2.1. The distribution of microlensed light-curve derivatives 

We are interested in the independent fluctuations in the 4 individual images of 
Q2237+0305. It is advantageous to look at the rates of change of the differences 
between image magnitudes because intrinsic fluctuation (time-delays are less 
than 1 day (eg. Schneider et al. 1988)) and the systematic component of error 
are removed. 

The basic tool of our analysis is the distribution of difference light curve 
derivatives (Wyithe, Webster &; Turner 1999a,b,c) in which the transverse ve- 
locity is a scaling factor. The cumulative histogram has a smaller value if the 
transverse velocity is higher. Thus by considering the shape of the cumulative 
distribution of light curve derivatives, information on the transverse velocity can 
be obtained. 

2.2. Monitoring data 

By combining the data of Kent & Falco (1988), Schneider et al. (1988), Irwin 
et al. (1989), Corrigan et al. (1991) and 0stensen et al. (1996), a light curve 
is produced for each of the four images. The most complete light curves are 
in R-band. There are a total of 61 published data points for each image taken 
between 1985, and 1996. To minimise the noise in the data, we averaged adjacent 
observations leaving 26 points. We simulate errors in the model light curves 
according to a Gaussian distribution. Am = 0.02 is taken as the 2a and la 
level in images A/B and C/D respectively. The level of noise in the sample is 
estimated by calculating the fluctuations that result from the application of the 
observational uncertainty and sampling rate to a flat light curve. Figure || shows 
the cumulative histogram of difference light-curve derivatives for Q2237+0305 
together with cumulative histograms for the mean, it la and ±2(7 cumulative 
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Figure 2. Left: Pv{Vgai < Vefj) and Right: U {{m) < M) under 
assumption of flat prior (dark line) and logarithmic prior (light line) . 

histograms of derivatives due to observational error only. The difference between 
these histograms demonstrates a statistically significant level of microlensing. 

2.3. Limits on effective transverse velocity 

We define a set of effective transverse velocities {14//}. 5000 simulations of 
the monitoring data were produced at each effective transverse velocity. An 
average histogram is also computed from each set of mock observations. For 
each mock observation at each pre-defined effective transverse velocity VJ=//, a 
mock measurement of effective transverse velocity fg// is made by minimising 
the KS difference between the average histogram for Veff and the histogram of 
the mock observation. Thus we calculate the function representing the likelihood 
Plhoodivef f\yef f) for observing fg// given an assumption for the true value (14//). 
Similarly, we find the effective transverse velocity (vobs) that best describes the 
observed microlensing rate. Using Bayes' theorem we calculate the posterior 
probability that the effective galactic transverse velocity Vgai is less than an 
assumed value K//: 

PviVgal <Veff\Vobs) = PlhoodiVobs\V^f f)Pprior{V^f f) clV^j f . (1) 



Our results are consistent for the assumptions of flat {pprior{Veff) oc dVeff) 
and logarithmic {pprioriVef f) oc ^^) prior probabihties for 14//. Pv{vgai < 
yeff\vobs) is shown in the left hand plot of figure The model presented here 
and in section |3] assumes a trajectory parallel to the C-D axis and no smooth 
matter. 14// is limited to be less than IOOOa/ (m) km sec~^ with 95% confidence. 

3. Determination of lower limit to microlens mass 

The velocity dispersion in the centre of the lensing galaxy is cr* ^ 215 km sec^^ 
(Foltz et al. 1992). This internal motion produces a microlensing rate oc (m)~2 , 
thus there is a minimum (m) which is consistent with observations. 
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3.1. Microlensing due to random proper motions 

We define the equivalent transverse velocity Vequiv as the effective transverse 
velocity of a static screen of point masses which produces a microlensing rate 
closest to that of the same model with Vequiv replaced by the combination of 
a random velocity dispersion (cr^,) plus a transverse velocity {Vtran) (Wyithe, 
Webster & Turner 1999a, b). The microlensing rates are considered equivalent 
when the histograms are most similar (have the minimum KS difference). Note 
that this process can be inverted to find the Vtran corresponding to a measured 
limit on K//. 

3.2. Limits on the microlens mass 

Vequiv quantifies the minimum level of microlensing that must be observed for an 
assumed physical transverse velocity of the lensing galaxy. Since the measured 
effective transverse velocity scales with (m) the cumulative probability for (m) 
is 

U l^m) <M= (^^"^ ) = / " ^ (^»"' < ^e//VM)) pp„„, {V,ff) dVtran 

(2) 

(Wyithe, Webster & Turner 1999c), where the priors PprioriVeff) are those used 
for the computation of PyiVgai < VeffVM)- Note that this distribution (shown 
in the right plot of figure y for the two assumed priors) does not represent the 
mass function, (m) is greater than O.IMq (95% confidence), and has a most 
likely value consistent with that of low-mass stars. 
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